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ABSTRACT 

Deep infra-red observations and long-term monitoring programs have provided dynamical evidence 
for a supermassive black hole of mass 3 x 10 6 M Q associated with the radio source Sagitarrius A* at 
the center of our Galaxy. The brightest stars orbiting within 0.1 parsecs of the black hole appear to 
be young, massive main sequence stars, in spite of an environment near the black hole that is hostile 
to star formation. We discuss mechanisms by which stars born outside the central parsec can sink 
towards the black hole and conclude that the drag coming from plausible stellar populations does not 
operate on the short timescales required by the stellar ages. We propose that these stars were dragged 
in by a second black hole of mass ~ 1O 3_4 M , which would be classified as an intermediate- mass 
black hole. We discuss the implications for the stellar populations and the kinematics in the Galactic 
center. Finally we note that continued astrometric monitoring of the central radio source offers the 
prospect for a direct detection of such objects. 

Subject headings: 
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1. INTRODUCTION 

Near-infrared imaging observations with speckle and 
adaptive optics techniques now allow the study of the 
central parsec of our Galaxy at an angular resolution 
of 0.05 arcseconds. Monitoring over the course of the 
last decade has provided the proper m otions for many of 
the stars orbiting in the central parsec ijEckart fc Genzefl 
Il997t lOhez et all 11998ft . Measurement of significant 
deviations from linear motion has yielded orbital so- 
lutions for several of the stars closest (< 0.016 pc) 
to the s upermassive bl ack hole (SBH) at the Galac- 
tic centerdEckart et al.l2002HGhez et al.ll2000l l2003albl : 
Schodcl ct al. 2002). While the older, fainter stars ex- 
hibit an apparently isotropic distribution of stellar veloc- 
ities, there exist separate populations of young stars, the 
motion of which appears to be dynamically unrelaxed. 
At distances ~ 0.1 parsecs from the black hole lie a group 
of He emission-line stars (the IRS 16 clust er) belonging 
to a t angcntially-anisotropic orbital family l|Genzel et alJ 
2000). These stars, with K band magnitudes ~9-14, are 
believed to be evolved, supergiant or Wolf-Rayet stars, 
with masses estimated at 30 — 100 Mq after correcting 
for the distance and the e xtinction towards the Galactic 
center l|Ghez et al.ll2003b|) . As such, these stars can only 
be ~ 1 — 10 Myr old. Closer to Sgr A* there is a coeval 
but apparently kine matic ally distinct popula tion of stars 
JGenzel et alJll99nl2003lGezari et alJl2002t lGhe"z et all 
l2003albl h This group is claimed to be radially anisotropic 
and of an earlier, O/B spectral type. It is the provenance 
of these two groups of stars and the origin of their pecu- 
liar kinematics that we wish to address. 

2. PREVIOUS WORK 

In situ formation of young massive stars is unlikely 
inside the central parsec as the tidal forces render it 
difficult for a sinking molecular cloud to survive long 
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enough to form stars close to the SBH. Clouds that 
are sufficiently dense to resist tidal shear, n > 10 8 
cm" 3 , are Jeans-unstable and fragment before they sink 



(Vollmcr & Duschl 2001). An alternative suggestion 
ijLevin & Bcloborodov 2003) is that the stars formed 
in an extended, self-gravitating gaseous disk that had 
in the past existed close to the SBH. However there is 
no direct evidence to corroborate this hypothesis. To 
explain the non-coplanar o rbits of the central cluster, 
ILevin fc Beloborodovl l|2003|) propose that the orbits of 
the closest stars have been affected by Lens-Thirring pre- 
cession due to an SBH spin that is not aligned with the 
disk axis. 

If the young stars were not born in their cur- 
rent positions, they must have migrated inward 
from larger Galactocentric radii. Massive star clus- 
ters <IGerhardll200lt IMcMillan fc Portegies Zwarll2003l 
Portegies Zwart. McMillan & Gerhard 2003) can sink 
via dynamical friction within Myrs but are tidally dis- 
rupted at a distance > 1 pc from the SBH where the 
specific gravitational binding energy is still ~ 100 times 
smaller than that of the most bound He-line star SO- 
2. Clusters impinging on the Galactic center on nearly- 
radial orbits are dispersed at their first passage near the 
SBH and result in a population of plunging, low-binding 
energy orbits unlike those of either the He line stars or 
of the central cluster. 

Stars that form in the molecular clouds of the circum- 
nuclear disk (CND) at ~ 1 — 2 pc, as well as those that 
are deposited by the clusters disrupted at ~ 1 pc, dif- 
fuse toward orbits of larger binding energy on the local 
relaxation time scale: 



T rel ~3 x 10 9 yr ( ^ 

U00kms-y 



( M* 



2 x 10 5 Af Q pc- 3 



\3M Q 

lnA x 
10 



-T7T , (1) 



where a is the local linear stellar velocity dispersion, M* 
is the average stellar mass, p is the local stellar density, 
and ln(A) is the Coulomb logarithm. Not only does this 
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exceed the age of the observed stars, but relaxation via 
star-star scattering would have also erased the observed 
kinematic peculiarities. 

The fundamental reason for the long relaxation time 
is that the stellar velocity dispersion a oc i? -1 / 2 diverges 
as the Galactocentric radius R decreases. Large veloc- 
ities weaken the effect of gravitational focusing in stel- 
lar scattering and decrease the average amount of en- 
ergy transferred in a single encounter. This limitation 
applies equally to all stellar mass objects. Only ob- 
jects significantly more massive than a star can sink on 
the timescales required. Star clusters, however, are not 
dense enough to survive intact in the strong tidal field 
l|Portegies Zwart. M cMillan fc Gerhard! 12003ft . Thus, 
the only astrophysical entity both massive and dense 
enough to satisfy the requirements is an intermediate- 
mass black hole (IBH). The orbital decay of an IBH 
is driven by dynamical friction on a shorter timescale 
Idf ~ (M*/Mi BH )T rc i - 1 - 10 Myr for an lO 3 - 4 Af 
IBH. Furthermore, it has been argued that such black 
holes form generically in dense, young stellar clusters as 
a result of the segreg ation of massive stars to the cluster 
center ()Spitzerlll969t) followed by the runaway merging 
in stellar collisions llPort egics Zwart & McMillan 120021: 
iRasio. Freitag fc Gurkanll2003jK 

3. STELLAR ORBITAL MIGRATION 

A small fraction of the parent stellar cluster may re- 
main bound to the IBH as it sinks in the Galactic po- 
tential, namely those stars originally located within the 
dynamical radius of influence of the IBH. Tidal stripping 
and ejection due to strong stellar encounters gradually 
remove stars from the cluster. In the absence of internal 
dynamical evolution, cluster stars orbiting at distance r 
from the IBH will be lost when they slip over the clus- 
ter Roche limit r = (Mi B H/AfsBH) 1/3 -R, where R is the 
distance to the SBH. 

However, stars are also ejected by scattering off 
other stars in t he cl uster. Using equation [37] of 
iLin fc Tremaind <ll98fl). the ejection timescale from a 
power-law cusp l|Bahcall fc WolJll976j) dominated by the 
IBH is 

\ 3 
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where tr c i, p c \ are the velocity dispersion and stellar mass 
density of the parent cluster at the original dynamical 
radius of influence of the IBH. 

The stars lost from the cluster are deposited over a 
range of radii. If stars are removed by tidal stripping 
alone, the profile of deposited stars reflects the original 
profile of the cluster, p(R) ~ i?~ 7 , where the Bahcall- 
Wolf value is 7 = 7/4. The stars most tightly bound to 
the IBH will be lost at a distance from the SBH of 
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Therefore to deposit the most bound star SO-2 at 
R ~ 0.01 pc, Mibh ~ 4 x 10 3 M Q is required. 



If the IBH orbital decay is ecce ntric and self-similar 
ijValtaoia. Valtonen fc Bvrdlll989| ). the profile may not 
be smooth; the young stars will be deposited in batches 
at a discrete set of locations corresponding to successive 
pericenter passages of the IBH. 

Even stars not directly bound to the IBH may be trans- 
ported if their orbits come close that of the IBH. The dy- 
namics of the interaction between stars and the SBH- IBH 
binary is very similar to the interaction of comets with 
the Sun- Jupiter binary. In the case of the Jupiter family 
comets, they are observed in the inner solar system after 
being scattered b y Jupiter from orbits with much la rger 
semi-major axes ijQuinn. Tremaine fc Duncanlll99flK In 
a similar fashion, successive weak encounters with the 
IBH cause a random walk in the orbital parameters, and 
this can nudge the peribothra of some field stars inward. 
One difference is that the IBH is gradually moving in- 
wards, in a fashion similar to planetary migration, for 
which planetesi mal/comet scattering may also be a con- 
tributing factor ijMurrav et alJ lT998) in some extrasolar 
systems. For this mechanism to be significant, the char- 
acteristic time /u _1/ ' 3 T or b between close encounters with 
the IBH must be shorter than the time yU 1//3 Tdf in which 
the IBH migrates a distance equivalent to the radius of 
its own sphere of influence, where T or b is the stellar or- 
bital period and p = A/ibh/A/sbh- This holds as long 
as 
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where we have assumed the stellar density p oc 1 / R 2 on 
these scales. Thus, the orbital decay of the IBH con- 
tinues to push some stars inward even after the origi- 
nal cluster has been tidally disrupted. The efficiency of 
these processes, however, is low, because most scattered 
stars are ultimately ejected — producing the Port cloud - 
in the case of comets l|Duncan. Quinn fc Tremainelll98'7l 
iFernandez fc Irjlll981(l — and a large population is needed 
at the outset. On the other hand, for an IBH born of a 
lO 5 Af cluster, the efficiency of the process need only be 
10~ 3 to explain the handful of IRS 16 stars. 

Dynamical friction drags the IBH toward the SBH 
until it reaches a critical separation where the bind- 
ing energy in stars with peribothra smaller than the 
SBH-IBH separation is the same as that of the IBH 
itself. Using the most recent determinatio n of cen- 
tral stellar cusp mass bv iGenzel et al.l (|2003), M cusp ~ 

1.3 x lO 4 M (i?/l") 163 , we infer this critical radius to be 
R ~ 0.2" - O.OO8(A/ibh/1O 3 M ) - 61 parsecs. It could 
be smaller if there is a significant dark stellar population 
interior to this orbit. Subsequent orbital decay proceeds 
at a decreased rate contingent on how efficiently star- 
star scattering and other forms of orbital diffusion feed 
stars into the emptied region (usually called the "loss 
cone" ) , thereby prov iding new material for slingshot ejec- 
tion by the binary llBegelman. B landford & Rccs Il980l 
IMilosavlievic fc Merritti 12002 lOuinlanl 119961) . The bi- 
nary shrinks by an octave for every MiBH-worth of stars 
that are removed. Stars are usually ejected by the cu- 
mulative effect of several encounters, with their orbital 
parameters undergoing a random walk. If the separa- 
tion decreases to 10 AU (AiiBH/lO 3 Af ), the emission of 
gravitational radiation dominates and will expedite the 
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coalescence of the black holes within a million years. 

4. THE DISTRIBUTION OF STARS IN THE INNER PARSEC 

If the He-line stars are dynamically associated with 
the IBH, then their orbits must either be similar to the 
orbit the IBH had at some, perhaps earlier, stage of 
its infall, or must at least still cross the latter if they 
are in the process of being ejected by repeated scatter- 
ing. To test this, we have examined the existing velocity 
data ijGenzel et alJ (2000); A. Ghez, private communica- 
tion) for Hc-line stars with three measured velocity com- 
ponents. Combined with the projected position, we have 
information on 5 of the 6 phase space coordinates neces- 
sary to compute the orbits, as well as complete solutions 
for a few close to the SBH(Ghc z et al.N2003a|) . To as- 
sess the importance of the unknown component (location 
along the line of sight, Z) we explored the range from 
Z = 0.5 R± to Z = 1.5 R±, where R± is the projected 
distance to Sgr A*. Figure ^ shows the resulting orbits 
which are plotted in a specific energy-angular momentum 
diagram. 




0.1 1 10 
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Fig. 1. — The filled circles show the Hc-line stars associated 
with the IRS 16 group. The error bars indicate the effect of 
varying the unknown z between z = 0.5r^ and z = 1.5r^. 
The open circles show the central cl uster stars for which com- 
plete orbits can now be determined I Ghez et al. 2003b) (thus 
these error bars are real). Triangles indicate approximate val- 
ues for central cluster stars for which complete orbits are not 
yet known. The solid line indicates the locus of orbits which 
cross an eccentric orbit with semi-major axis 1" and eccentric- 
ity 0.82. An IBH on this orbit can interact with all the stars 
shown. The dashed lines indicate the same for circular orbits 
at 2" and 0.2" instead. The dotted lines indicate the locus of 
circular orbits (rightmost) and orbits with eccentricity 0.9. Eq 
and Lq are calculated assuming an SBH mass of 3 X 10 6 Mq 
and a semi-major axis of 1" = 0.04 parsecs. 

The IRS 16 stars all cross a circular orbit located 2 
from Sgr A* . As such they are consistent with being the 
remnant of a disrupted IBH host cluster or, alternatively, 
are currently undergoing the scattering-induced diffusion 
and the eventual ejection if the IBH is located at this ra- 
dius. The more tightly bound stars, the Sgr A* cluster, 
do not cross this orbit. However, it is possible for an ec- 
centric IBH orbit to interact with both the IRS 16 stars 



and the Sgr A* cluster. The solid curve in Figure de- 
lineates that region of parameter space (that part below 
the curve) in which stars interact with an orbit of semi- 
major axis of 1 (0.04 pc) and eccentricity 0.82, which 
covers all the early-type stars. 

It is also possible that the IBH orbit is indeed not very 
eccentric, and that the central cluster stars once did cross 
the IBH orbit, but have since then received a significant 
perturbation near peribothron that resulted in a new ap- 
bothron within the orbit of the IBH. The star is then 
trapped until further IBH inspiral brings the orbits into 
conta ct again. The claim ed radial anisotropy of these 
stars llSchodel et al.ll2003|) is consistent with such a sce- 
nario. The perturbation could have come from a grazing 
stellar collision or from the dynamica l disruption of a 
binary system ijGould fc Quillenll2003|) . For the pertur- 
bation to be strong enough to trap the star in a closer-in 
orbit, the collision must be grazing because the veloci- 
ties in this region are comparable to the escape velocities 
from the stars themselves, i.e., the dynamical evolution 
is dominated by physical collisions rather than two-body 
relaxation. 

We find that the spectral differences between the cen- 
tral cluster stars and the IRS 16 stars lend support to 
this picture. Although both are young and massive, the 
former appear to be main sequence stars of roughly O 
or B type, while the latter are extended, evolved stars 
in the supergiant or the Wolf-Rayet stage. The more 
compact nature of the central cluster stars may indicate 
that the outer envelopes have been stripped in grazing 
collisions. Not e that this is d ifferent from the alternative 
proposal (e.g-. lGenzel et al]l2003j) in which the observed 
stars are assembled from lower mass objec ts in stellar 
merger s. The normal rotation rate of SO-2 ijGhez et alJ 
2003b) and the potential radial anisotropy of the orbits 
are both consistent with a limited role of collisional in- 
teractions. 

If an IBH is present at the Galactic center, there is also 
the possibility of directly observing high velocity stars 
that have been ejected from the region. This would es- 
pecially be true if the SBH-IBH binary has hardened past 
the stalling radius and is continuing to evolve slowly via 
the gravitational slightshot of stars that diffuse into the 
loss cone. The velocities of the ejected stars will reflect 
the specific binding energy of the IBH at the time of 
ejection. 

Direct constraint on the IBH hypothesis may be pro- 
vided by the non-inertial motion of Sgr A*. The pres- 
ence of an IBH at the Galactic center will be revealed 
in the apparent proper motion of the radio image of 
Sgr A* relative to the mean ~ 6 mas y r -1 solar drift 
{Backer fc Srameklll999t iReid et al.lll999j) . Assuming a 
circular orbit of the IBH, the amplitude of the grav- 
itational reflex of SBH relative to the barycenter of 
the system as a whole, ^ref ~ Mibh/^sbh where 
\& is the angular separation between the black holes, 
must be larger than the cumulative positional resolution 
A 1 !* < 1 mas of the radio telescope to achieve detection. 
In turn, the total projected distance traversed by the 
SBH, equal to * rcf ~ T^ImuR^iG/ RM SB n) 1/2 for a 
monitoring program of duration T m and distance to the 
Galactic center R gc ~ 8 kpc, must also be larger than 
A^J*. These constraints are summarized in Figure |21 
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5. CONCLUSIONS 

In this paper we have proposed a model that ad- 
dresses the peculiar nature of the young stars in the 
Galactic center. Most importantly, the presence of an 
IBH can deliver the massive young stars to their ob- 
served location within the timescale required. In ad- 
dition, our scenario provides a way to link the two 
disparate groups of stars (the IRS 16 group and the 
Sgr A* cluster) within a single evolutionary scenario. 




M IBH ( M ©) 

Fig. 2. — Assuming a circular orbit around an SBH of 
3 X 10 6 Mq, we can rule out an IBH with mass Mibh and semi- 
major axis a by measurement of an astrometric wobble of the 
radio image of Sgr A* . The shaded regions show the detection 
thresholds for astrometric resolutions of 0.01, 0.1 and 1 milliarc- 
seconds, respectively, assuming a monitoring period of 10 years. 
The dashed lines indicate coalescence due to gravitational radi- 
ation in 10 6 and 10 7 years, respectively. 

In this scenario, once the IBH inspiral stalls, the stellar 
orbits continue to evolve under the action of scattering 
during close encounters with the IBH. Most stars will 
eventually be ejected, but some stars may be trapped 
near the SBH if they receive a sufficiently large perturba- 
tion (such as due to a grazing stellar collision). This may 
also serve to generate the observed spectral differences 
between the extended and apparently undisturbed He- 
line stars, and the more compact but apparently contem- 
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poraneous Sgr A* cluster stars. Furthermore, the tangen- 
tial anisotropy of the He line stars may reflect their still 
somewhat limited dynamical evolution under scattering, 
whereas the claimed radial anisotropy of the Sgr A* clus- 
ter stars may reflect their having diffused to highly eccen- 
tric orbits before being scattered and trapped. However, 
this radial anisotropy is still uncertain because its deter- 
mination is subject to severe selection effects. 

Some questions remain open: Although the mass den- 
sity of luminous stars can be measured, the true den- 
sity profile of the central parsec is not well-known. Our 
adoption of the observed density profile is conserva- 
tive as there are good reasons to believe that IOMq 
black holes will have segregated within few relaxation 
times to dominate the density of the cent ral cluster 
ijMiralda-Escude fc Gouldll200a IMorrisl ir993L If such a 
dark cusp exists the central density could be as large as 
lO 9 M pc -3 and thus the stalling radius of the IBH could 
be much smaller than estimated above. Alternatively, if 
the IBH infall is an episodically recurring phenomenon, 
then the injection of a steady stream of IBH into the cen- 
tral parsec would help to maintain an evacuated region 
in the stellar distribution. 

Finally, IBH are implicated in the form ation of SBH in 
galaxies in a vari ety of wavs ( | R eesj|1984fl . Indeed, it has 
been suggested (Ebisuzaki et al.ll2001|) that SBH form 
from the collapse of a cluster of IBH. It is intriguing that 
to gradually accumulate the present SBH at the Galactic 
center by this process requires IBH captures at the rate 

r ~ 3xl ° y1 ' Ux106mJ U^oJ ' (5) 

With this rate one would expect the most recent entrant 
to be a few Myr old, which coincides with the time scale 
derived above from the stellar ages. Such a process could 
have a significant influence on t he history of galactic nu- 
clei lIHughes fc Blandfordll2fK)^ . 
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